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ABSTRACT 

We present a path-integral likelihood formalism that extends parameterized likelihood 
analyses to include continuous functions. The method finds the maximum likelihood 
point in function-space, and marginalizes over all possible functions, under the assump- 
tion of a Gaussian-distributed function-space. We apply our method to the problem 
of removing unknown systematic functions in two topical problems for dark energy 
research : scale-dependent galaxy bias in redshift surveys; and galaxy intrinsic align- 
ments in cosmic shear surveys. We find that scale-dependent galaxy bias will degrade 
information on cosmological parameters unless the fractional variance in the bias func- 
tion is known to 10%. Measuring and removing intrinsic alignments from cosmic shear 
surveys with a flat-prior can reduce the dark energy Figure-of- Merit by 20%, however 
provided that the scale and redshift-dependence is known to better than 10% with a 
Gaussian-prior, the dark energy Figure-of-Merit can be enhanced by a factor of two 
with no extra assumptions. 

Key words: Cosmology: theory - large-scale structure of Universe 
Methods: data analysis - analytical - statistical 



1 INTRODUCTION 

The Standard Cosmological Model is a phenomenological 
model containing two unknown components, dark energy 
and dark matter, and two untested assumptions, general 
relativity on cosmic scales and cosmological inflation. To 
understand this model at a deeper level, by testing the evo- 
lution of dark energy, the nature of the dark matter, sig- 
natures of inflation and possible deviations from Einstein 
gravity, we need highly sensitive cosmological probes that 
must constrain an expanded set of cosmological parameters 
to (sub-)percent accuracy. In such a scenario the correct as- 
sessment of systematic effects is of critical importance. 

To account for systematic effects it is common to use a 
parameterized description of the effect where each nuisance 
parameter is calibrated to some accuracy and removed. How- 
ever, there are many cases where characterization of the 
systematic effect by a parameterized function is not well- 
justified, and in some cases the systematic is effectively an 
unknown function. The assumption of a given functional 
form can lead to a bias in measured cosmological param- 
eters, or at least an underestimate of the contribution to 
errors due to uncertainty in the systematic function. 

In Taylor & Kitching (2010) we described a new method 
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which allowed the analysis of large-dimensional parameter- 
spaces composed of cosmological parameters and nuisance 
parameters describing residual systematics, in both cal- 
ibrated (Gaussian prior) and self-calibrating (flat-prior) 
regimes. Here we present a new method that can account for 
systematics in a model-free way such that the impact of ev- 
ery possible function on the cosmological signal is accounted 
for. To do this we introduce a path-integral marginalization 
approach in which we maximize the likelihood and marginal- 
ize over a space of functions rather than parameters. We 
construct a likelihood functional in parameter space that 
accounts for the impact of all known systematic functions. 
The function-space may be either weighted with a flat prior 
about a flducial function, where the data itself is used to 
measure and remove the systematic effect, or if there is ex- 
ternal data weighted with a Gaussian prior, with a variance 
that is zero when the systematic function is known. 

We derive an expression for the Fisher matrix that ac- 
counts for functional marginalization, and apply this method 
to two test cases in cosmology in which the cosmological sig- 
nal is contaminated by a systematic effect with a poorly un- 
derstood functional behaviour. The flrst is scale-dependent 
galaxy bias (e.g.. Peacock & Smith, 2000; Conway et al., 
2005; Hamann et al., 2008; CressweU & Percival, 2009) in 
redshift surveys, an unknown function that relates the un- 
derlying dark matter distribution to the galaxy distribution. 
The second is cosmic shear intrinsic alignments (e.g. Brown 
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et al., 2002 ; Crittenden et al , 2001, Catelan et al„ 2001; 
Heymans et al., 2004; Hirata & Seljak, 2004; Kitching et 
al., 2007; Bridle & King, 2007; Joachimi & Bridle, 2009; 
Joachimi & Schneider, 2008, 2009; Mandelbaum et al., 2009; 
Schneider & Bridle, 2009; Kirk et al., 2010) in weak lensing 
surveys, a contaminant with a complex astrophysical func- 
tion that relates the local density to the ellipticity of back- 
ground galaxy images. In both cases we show how cosmo- 
logical parameter constraints are expected to degrade when 
we marginalize over all systematic functions. 

This work is an extension of the methods of Taylor & 
Kitching (2010) in which a formalism for analytic marginal- 
ization was developed for a finite number of parameters, and 
is related to the Functional Form-Filling method of Kitching 
et al. (2009) in which a numerical approach to assessing the 
effects of unknown functions was presented. 

This paper is organized as follows, in Section[2]we intro- 
duce path integral marginalization, in Sections and U we 
apply the method to galaxy bias and cosmic shear intrinsic 
alignments respectively. In Section[5]we present conclusions. 
Throughout we have put technical detail in a series of Ap- 
pendices. 



2 PATH-INTEGRAL MARGINALIZATION 

We consider an arbitrary likelihood functional, L[G, tp{x\6)], 
which is a function of a finite number of cosmological pa- 
rameters of interest, 6, and a finite number of continuous 
nuisance functions, xpa{x\G), labelled by Greek subscripts. 
The nuisance functions could depend on position x — x, ot 
wavevector a; = fc if we are working in Fourier space. The 
nuisance functions may themselves also depend on the cos- 
mological parameters. For convenience we shall work with 
the log-likelihood functional; 



c[e,%p{x\e)] = -2inL[e,%p{x\e) 



(1) 



Expanding the log-likelihood functional to second order 
around an arbitrary nuisance function we find 
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(3) 



where 5%l)a{x\0) is the displacement in function-space from 
the fiducial function and the maximum-likelihood solution, 
and we have assumed the inverse of the second functional 
derivative of the log-likelihood exists. 



2.1 Flat Priors 

Marginalizing over all paths for the functions, assuming a 
fiat-prior, with boundaries sufficiently far away they do not 
contribute to the integral, we find 



£ = £(,+ I dx'dx'Ti In 
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where we have omitted an unimportant constant. 

The first and third terms in this expression represent 
the maximum log-likelihood value, which can be seen by 
substituting the solution for the nuisance function tl>{x\9) 
given by equation ^ into equation 0. As we always find 
the maximum-likelihood solution for the nuisance functions, 
the marginalized likelihood is independent of the fiducial 
nuisance functions, if}{x\9) in the Gaussian approximation. 
The second term in equation Q is the width of the likeli- 
hood function. 

In Taylor & Kitching (2010) we showed that in the 
case of marginalization over discrete nuisance parameters 
the Fisher information is preserved about the cosmological 
parameters of interest. This result holds in the case of con- 
tinuous functions, '0(x|0). 

2.2 Gaussian Priors 

If the nuisance functions are constrained by another, in- 
dependent experiment we can include this constraint as a 
Gaussian prior on ■tp{x\0), with covariance matrix 



Cc.f,{x,x') = {s^p^{x\e)Si>0{x'\0)). 



(5) 



Including this constraint in the marginalization path integral 
we find 



C = Co - - / dx'dx' 

2 / Stpa{x 



-M^i^ix ,x )- 
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Mal3{x',x") 



S'Co 



+ 2C-',{x',x") 



(7) 



6^a{x')5llji3{x" 

and we assume the matrix Mo-a is invertible. The fiat-prior, 
which we have assumed is non-zero over a finite but suffi- 
ciently large region, is not a limiting case of the Gaussian- 
prior. When the covariance matrix of the Gaussian-prior 
goes to infinity, the prior distribution becomes arbitrary 
close to zero everywhere, and the infinite wings of the prior 
suppress information. 



2.3 Gaussian Likelihoods with Systematics in the 
Mean 

A common distribution for cosmological datasets is that of a 
multivariate Gaussian where, as discussed in Taylor & Kitch- 
ing (2010), the effects of cosmological parameters and nui- 
sance parameters can appear in both the mean and/or the 
covariance. The likelihood Gaussian function is given by 

Co = ADC'^AD' + Tr In C, (8) 

where AD — D ~ fi is the variation of the data around 
its mean value, /j. — fi[d,xl>{x\0)], and C[9,ip{x\0)] — 
(ADAD^) is the data covariance matrix. 

For simplicity we shall assume the nuisance functions 
and cosmological parameters are only in the mean term, fj,, 
however as shown in Taylor & Kitching (2010) results can be 
generalized for nuisance functions in the covariance matrix. 
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We assume the curvature of the hkelihood in function-space 
can be approximated by its ensemble average, so that 

1 / 5'^Co{x,x') 



Fal3{x,x') 



= Tr C 



2 \5i}6{x)5iPfi{x') 

^5^J,['^{x)]S^l^[^|^{x)] 



(9) 



5i>a{x) Sil;is{x') 
is the generalized Fisher matrix, and 

Si,ix') ^^'~ 

is the gradient of the log-hkelihood. The maximum- 
likehhood solution for the function is 



5ip^{x) = ADC^ / dx' 



Sip 13 (x' 

e log 
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where again Greek indices represent a discrete infinite set of 
functions. 



2.3.1 Flat prior 

We can now write the full path-integral marginalized log- 
likelihood for a flat-prior as 

Co = ADC M AD' + Tr In Cm, (12) 

where 

Cm = [C-'' -C-'^PC-'T'' (13) 

is the new marginalized covariance matrix, 



(14) 



and we assume the functional Fisher matrix, Fa/3{x,x'), 
is invertible. The marginalized covariance matrix. Cm = 
Cm (9), is now a function of cosmological parameters 
through the inclusion of the functional derivatives of the 
mean, ft{6). 

2.3.2 Gaussian prior 

In the presence of a Gaussian-prior, the form of the 
marginalized likelihood and covariance matrix still holds, 
with the matrix P now given by 

Here the Map is a sum of the path-integral Fisher and a 
prior covariance term 

Maf}{x,x') = Faf}{x,x) + C~l{x,x). (16) 

This contains information on the nuisance parameters from 
the prior Cai3{x,x') and the data itself Fai3{x,x'). If we in- 
clude this prior covariance on the nuisance functions we find 
(see Appendix A for details) that the marginalized likeli- 
hood is still a Gaussian in the data, with a new covariance 
matrix; 

This form of the marginalized data covariance matrix has 
the clear advantage that it does not require the inversion of 
a functional Fisher matrix, or prior covariance matrix. 



2.4 Marginalized Fisher Matrices 

The Fisher Matrix for the marginalized likelihood is given 
by (Tegmark, Taylor & Heavens, 1997) 

Fab = ^Tr K/i^ + MaMblCM , (18) 

where Roman indices, {a,b), denote cosmological parame- 
ters, and for path-integral marginalized covariance Cm is 
given by equation ((13} for a flat prior and equation (|17|l for 
a Gaussian prior. We have implicitly taken account of the 
cosmological parameter information contained in the mean 
of the data via fJ,{d). However, even though the marginal- 
ized covariance matrix. Cm, now depends on cosmological 
parameter through the functional derivatives of the mean, 
fi, at the level of the Fisher matrix Cm is a constant in 
parameter-space. 

Throughout Section 2 we have made the assumption of 
a Gaussian-distributed function-space. If this assumption is 
dropped then the functional integration can be done numer- 
ically, one way to do this integration is using Monte-Carlo 
integration techniques, or functional form-filling (Kitching 
et al., 2009). We discuss these alternatives in Appendix B. 

Having developed the formal aspects of path-integral 
marginalization to cosmological likelihoods and Fisher ma- 
trices, we now apply it to the case of marginalization over 
two cases of systematic effect in Cosmology. The first is 
scale-dependent galaxy bias, and the second is the case of 
intrinsic alignments in weak gravitational lensing. 



3 SCALE-DEPENDENT GALAXY BIAS 

As an example let us consider the case when the mean signal 
is an estimate of the galaxy power spectrum measured from 
a galaxy redshift survey; 



fi{k,r) = P,,{k) + N{r), 



(19) 



where N{r) = l/n(r) is the distance-dependent galaxy shot- 
noise and n{r) the mean galaxy number density. The galaxy 
power spectrum is a biased representation of the underlying 
matter power spectrum 



Pgg{k)=h\k)Pss{k), 



(20) 



where Pss{k;) is the matter power spectrum and b{k) is a 
dimensionless, scale-dependent galaxy bias function. This 
scale-dependent galaxy bias only depends on scale and not 
phase, so is a function of fc = |fe| only. While the bias func- 
tion is, in this model, deterministic, it is in general com- 
pletely degenerate with the matter power spectrum. Hence 
we need to introduce a Gaussian prior from the start, where 
the prior knowledge could come from weak lensing mea- 
surements (e.g. van Waerbeke, 1998; Bernstein, 2009) or a 
semi-analytic galaxy bias model. In either case there will be 
some level of uncertainty associated with the bias model, 
5tp{k) — 5b{k). We can, of course, extend this model to in- 
clude stochastic bias effects (e.g., Dekel & Lahav, 1999). 

In the following we include the effect of linear redshift- 
space distortions by the transformation (Kaiser, 1987) 



P^Ak) 



[b{k) + ff^lfPss{k), 



(21) 



where / = d\n5,n/d\na is the growth rate of matter per- 
turbations, and ^k is the cosine of the angle between the 
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wavevector k and the line-of-sight. We shall assume we are 
working on scales large enough to ignore non-linear finger- 
of-god effects. If we assume the power is a Gaussian random 
variable, the covariance of the galaxy power is given by 

C, = 2[P^,ik) + N{r)]\ (22) 

Applying equation ([17]), assuming a Gaussian prior, the 
functionally marginalized covariance is 



CM{k) = 2[P;,(fc) + N{r)Y + B{k), 
where 



Bik) 



Pbjk) 
b2{k) 



\Pooik)\ 



(23) 



(24) 



is the increase in the covariance matrix after we have 
marginalized over all bias functions. As the galaxy bias is 
independent of phase, and both the galaxy and matter dis- 
tributions are statistically homogeneous, we only need to 
know the power, or functional variance, in the bias; 



Ptik) = {\Sb{k)\'), 



(25) 



were the averaging is taken over all functions of b{k). 

If we substitute the redshifted galaxy power spectrum 
into the Fisher Matrix formalism, (see e.g., Taylor & Watts, 
2001; Seo & Eisenstein, 2003; Burkey & Taylor, 2004), we 
find the Fisher Matrix including path-integral marginaliza- 
tion over all possible bias functions is 

d^fc din P^g{k) din P^g{k)^ 



Fab - 2 

where 



Kff(fe) 



(27r)3 



dOa 



-Kff(fc), 



[P^Ak)] 



2[P|,(fe) + iV(r)12 + B(fc) 



(26) 



(27) 



is the effective volume of the survey. As B{k) also depends 
on Pgg{k), the effective volume fully specifies our survey. In 
the limit of no systematic uncertainty in the bias, Pb{k) = 0, 
or B{k) = 0, this reduces to the usual Fisher matrix for the 
redshifted power spectrum. 

We demonstrate the effect of path-integration over all 
possible galaxy scale-dependent bias functions with a Gaus- 
sian prior in Figure [T] where we consider a Euclid spectro- 
scopic survey (Laureijs et al., 2009) with a survey volume of 
19.7 h~''Gpc'' and a maximum wavenumber of 0.5/iMpc~^. 
We shall assume that the signal-to-noise ratio of the sur- 
vey is unity, nPgg — 1 (Seo & Eisenstein, 2003). We use 
a cosmological parameter set that allows for curved cos- 
mologies with parameters Qrn, ^de, ^b, h, erg, and Us given 
by (0.25,0.7,0.75,0.8,0.95), and parameterize the dark en- 
ergy equation-of-state using a first-order Taylor expansion, 
w{z) — Wo + {1 — a)wa (Linder, 2003; Chevallier and Po- 
larski, 2001) with {wo,Wa) — (—0.95,0.0). We assume a 
fiducial bias function of b{k) = 1. In Figure [T] we show the 
marginal errors on Qm , wo and Wa respectively, as a function 
of the functional variance in the bias. In the limit that the 
bias variance is small the accuracy converges towards the 
systematic-free case. However, as the bias variance increases 
beyond Pi,(fc)/b^(fc) « 10~^, the uncertainty on cosmological 
parameters grows large. Using this result and from equation 
(|27p we define an approximate scaling relation 
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Figure 1. The marginalized errors on Qm, (lower, black line), uiq 
(middle, red line) and «)a (upper, blue line) for a Euclid spectro- 
scopic survey (Laureijs et al., 2009) marginalized over all possible 
bias functions b(k). We varie the functional variance Pi,(k)/b^{k) 
which weights the space of functions that are marginalized over 
about a fiducial function. In the limit that Pi,{k)/b^{k) — > oo 
all possible functions are marginalized over, in the limit that 
Pl,{k)/b'^(k) — ^ the bias is fixed at the fiducial function. 



for cosmological parameter constraints to be unaffected by 
uncertainty in the scale-dependent bias. 

Constraints on the bias function of Pf,(fc)/b^(fc) ~ 10~^ 
can be achieved using large area surveys and techniques such 
as shear-galaxy cross correlations (e.g. van Waerbeke, 1998; 
Bernstein, 2009) or bispectrum statistics (e.g. Matarrese et 
al., 1997; Verde et al., 1998; Guo & Jing, 2009). A more 
complete analysis should include the redshift-dependence of 
galaxy bias, as well as scale-dependence, and also consider 
stochastic bias effects. However, consideration of equation 
(27) suggests that all of these effects will only become im- 
portant if they become larger than the sampling and shot 
noise terms. 



4 COSMIC SHEAR INTRINSIC ALIGNMENTS 

We now consider the effect of intrinsic alignments (lA) on 
weak lensing. Cosmic shear uses the weak lensing distor- 
tion caused by large-scale structure along the line of sight 
as a probe of cosmology; this assumes that galaxies are 
randomly aligned, so that the mean observed ellipticity is 
zero. However there are two effects which act to align galax- 
ies by adding 'intrinsic' ellipticity which is not part of the 
cosmic shear signal. Intrinsic-Intrinsic (II) alignment is the 
small-scale effect of galaxies aligning due, for example, to 
local tidal forces. We will ignore this effect in our example 
since this can be removed by neglecting close galaxy pairs 
in angle and redshift from the cosmic shear signal (e.g Hey- 
mans & Heavens, 2003). The Galaxy, or shear. Intrinsic (GI) 
alignment (Hirata & Seljak, 2004) is the correlation between 
the foreground galaxy ellipticity and the induced shear in a 
background galaxy and this is more difficult to correct for. 
The tomographic cosmic shear power spectra C^ {£), 
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where the two-point correlation of the shear field is binned 
in redshift, can be written as 



C?^{1) 



GG/ 



drW:- {r)Pi 



I 



\Sk{r) 

where lensing weight can be expressed as 



/GG 



w^r (0 = 



Slir) ' 



and the kernel is 



qi{r) = 



2a(r) 



dr' pi{r 






(29) 



(30) 



(31) 



We follow the notation of Joachimi & Bridle (2009). The 
comoving distance is r, rn is the horizon distance, while 
Skir) — sin(r),r, sinh(r) for curvatures k = — 1,0,-|-1, a 
is the scale factor and Pss{k;r) is the 3D density-density 
matter power spectrum. The comoving galaxy probability 
distribution is given by Pi (r) . The ij subscripts refer to red- 
shift bins, where the shear field is approximated as a series 
of correlated 2D planes. This can be generalized to the full 
3D shear field (Kitching & Heavens, 2010). 

In a similar way the GI power spectrum can be written 
as an integral over the density-shear power spectrum; 



Cr/i£)= / drW,']'{r)Ps^ 



£ 



(32) 



where we have defined the GI weight function as 

W,f{r) = n,{r)q,{r)/Sl{r), (33) 

and ni{z) is the number density of galaxies in the redshift 
bin i. We can now write the mean signal as 

M. {£) = C%° {£) + C^',,yM + iV.j , (34) 

where Greek subscripts v — {i,j) denote pairs of redshift 
bins, and the notation j > i indicates that the matrix is 
non-zero when the lensed galaxy is behind the intrinsically 
aligned galaxy. The shear shot noise term Nij = 5ija1/2ni 
(Hu, 1999), where a^ is the variance in the intrinsic ellipticity 
distribution, only adds to the Cij term. The covariance of 
the mean is 

c^^ii) = cf^{i) [cfi°(e) + cZ>M + cZ<M] 

+ [C?°{1) + Cg(,,>m(^)] [cfi'^W + Cfi'iy^H)] . (35) 

Currently there are two approaches to removing the 
GI contaminant; nulling and modelling. In nulling the GI 
contaminant is removed from the data by transforming the 
data vector, by ordering the transformed modes in order of 
contamination, and removing the most contaminated ones 
(Joachimi & Schnieder, 2008, 2009). However, nulling also 
removes cosmological information and can bias the recovered 
cosmological parameters. 

The second approach is modelling, which uses a phe- 
nomenological (e.g. Heymans & Heavens, 2003) or theoret- 
ically motivated model (e.g. Schneider & Bridle, 2009) and 
marginalizes over the parameters in that model. The prob- 
lem with the modelling approach is that it is sensitive to 
the choice of model and the resulting cosmological param- 
eter errors are dependent on this choice (Bridle & King, 
2007; Kitching et al., 2008; Joachimi & Bridle, 2009). Fur- 
thermore numerically marginalizing over lOO's of parame- 
ters, required to characterize the model, is computationally 



expensive, although semi-analytic marginalization could be 
used to reduce this expense (Taylor & Kitching, 2010). 

An alternative, third way, which we propose here is to 
marginalize over all possible intrinsic alignment functions 
using path-integral marginalization. 



4.1 Intrinsic Alignment Removal 

In this Section we investigate two distinct cases. The first 
case is marginalization over the uncertainty in the density- 
ellipticity relation itself. In this first case we assume both a 
fiat-prior (the self-calibration case), and a Gaussian-prior. 
A Gaussian prior may use information on the density- 
ellipticity relation gained from, for example, galaxy-galaxy 
lensing. In the second case we make some assumptions about 
how the density-ellipticity power spectrum is related to the 
matter power spectrum via a set of intrinsic-alignment bias 
functions, and assume only a Gaussian prior. 



4-1.1 Intrinsic Alignment removal: flat prior 

The functional variation of the lA signal with the density- 
ellipticity cross-spectrum is given by (see Appendix C) 



5C?/.i>A^) _^^^Gi ,.2^5d{£-1') 



SP-,s{e'/r,r) 



w^^,.>.(0- 



£ 



(36) 



For a fiat-prior, where we use the data itself to fit and 
marginalize over the lA term, we use equations 1)13^ and 
(11411 to find the marginalized covariance 

C^iW = [C-J{£) - C'l, (^)P^v' (^)C:,t(^)]-i (37) 

where 



Pm.(^)= / drdr'Wl:\r)[F{£-r,r')]-^Wl;\r') 



(38) 



and the functional Fisher matrix for the lA power spectrum 



F{£-ry) = Y,C-l{£)W°'{r)W^'{r') 



(39) 



To invert the functional Fisher matrix we bin the radial 
GI weights in equations (|39p and l|37p using 200 radial bins, 
while the angular terms are always diagonal. The Fisher ma- 
trix for cosmological parameters is given by equation (|18p . 
this is more explicitly shown for the lA case in Appendix C. 

To make a qualitative assessment of the impact of in- 
trinsic alignment removal we consider a Euclid weak lens- 
ing survey (Refregier et al., 2010) of 20,000 square de- 
gree, with median redshift of Zm = 1.0 and 35 galaxies per 
square arcminute with a photometric redshift uncertainty of 
0.03(l-|-z). We assume the galaxy number density is given by 
n{z) (X z'^eyi]i{ — lAz/z,nY'^, and use the same cosmological 
model as in Section 3. 

To quantify the effect of path-integral lA removal on 
dark energy surveys we use the dark energy Figure-of- 
Merit (FoM; Albrecht et al., 2006), defined as the area con- 
strained in the uncorrelated {wo^Wa) plane with FoM = 

'^l[\f{F^LoFwLa - {F^iiLaf)]- We assume the fiducial 
model for the GI Intrinsic Alignment is given by the lin- 
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ear/nonlinear alignment model of Bridle & King (2007; see 
also Section liTJJ Fl. 

For comparison we find that, for a default model with 
fixed lA and no cosmological dependence of the lA term (i.e., 
cosmic shear only, treating fA as additional known noise 
term), the FoM = 130. If we allow the ellipticity-density 
relation to be fixed by the data, with a flat-prior and allow 
a cosmological dependence in the weak-lensing weight of the 
GI term, we find the FoM is degraded to FoM « 100 (see 
Figure [2} . In this case, information is gained because of the 
extra cosmological dependence of the GI weights, but also 
lost as the data now has to measure the density-ellipticity 
relation. 

It is interesting to compare the path-integral lA removal 
approach to the nulling method of Joachimi & Schneider 
(2008, 2009) which uses the radial information in the lens- 
ing kernel to null the GI term from the shear data vector. 
Nulling makes a number of assumptions about the lens and 
source planes being thin sheets, and fixes the cosmology 
in the lensing kernel. An effect of nulling is to introduce 
a bias in cosmological parameters due to its effect on the 
shear power. Nulling finds a factor ~ 2 reduction in a global 
FoM, combining errors from all cosmological parameters, 
which seems comparable to path-integral marginalization. 
We leave a full comparison of methods for future work. 



4-1.2 Intrinsic Alignment removal: Gaussian prior 

If, instead of self-calibration, we use external data to con- 
strain lA with a Gaussian-prior, we find the path-integral 
marginalized covariance, equation (|17|l . can be written as 
(see Appendix C) 



1000 



ctiW 



c^Ae) + 



dral{l/r,r)W^/{r)wSr^(r), 



(40) 



where a%{£/r,r) = {\5Pjs{£/r,r)\^) is the prior uncertainty 
in the shear-density relation. We again make no assumption 
about the relationship between matter and ellipticity. The 
path-integral Fisher matrix, marginalizing over all possible 
shear-density cross-power, is again given by equation (|18p 
and in Appendix C. 

In Figure [2] we show how the dark energy FoM changes 
as the fractional scatter in the ellipticity-density power spec- 
trum; a[P-,s{£/r,r)]/P.ys{£/r,r). The lower solid line extends 
our default model from Section 14.1.11 by treating the lA 
term as additional noise where the GI term is independent 
of cosmological parameters but the ellipticity-shear relation 
is allowed to vary (constrained by the prior on P.ys{£/r,r)). 
In this case the FoM is rapidly damped from FoM — 130, 
as the scatter is increased above a{P^s)/P'yS ~ 10""^. As an- 
ticipated in Section 2.2, as the scatter is allowed to increase 
this adds greater uncertainty to the model, with ever larger 
fluctuations in the ellipticity-density relation. As a result the 
information is actually damped below that of the flat-prior. 

The upper curve in Figure [5] allows the redshift- 
dependence of the GI weights to vary with cosmology and 
act as a source of information. This increases the FoM, in 
the limit of no uncertainty in the lA term, to FoM = 300, 

^ Our code is an extension of the iCosmo package (Refregier et 
al., 2008) and is available on request 
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Figure 2. The dark energy Figure of Merit (FoM) for a Euclid- 
like tomographic weak lensing survey (Refregier et al., 2010) 
marginalized over all possible non-linear alignment bias func- 
tions. We show results as a function of the fractional scatter in 
the Pg^iik) function-space, assuming a Gaussian prior. The lower 
black line (GG) shows the FoM from cosmic shear alone, assum- 
ing the GI term does not depend on cosmological parameters. 
The upper red line (GG+GI) shows the FoM when we also allow 
the GI rcdshift dependence to vary with cosmological parameters. 
The horizontal (blue) line shows the FoM for the self-calibration 
case in which flat prior in function-space is assumed (equation 
I37| l — in this case the shear and GI power spectra are dependent 
on cosmology (similar to the upper red line). 



showing that the known lensing dependence of the GI term 
contains useful cosmological information. Even though the 
FoM drops off rapidly, as the scatter is increased, it still 
remains a factor of a few greater than ignoring the cosmo- 
logical information in the GI term. 



4-1.3 The Intrinsic Alignment Model 

In this Section we will assume that the shear-density power 
spectrum, P-^s, is related to the matter power spectrum by 
some unknown ellipticity bias which is a function of redshift 
and scale (e.g., Bernstein, 2009) 



P-,s{e./r,r) = bi{£)bi{r)Pss{£/r,r) 



(41) 



and that the dimensionless bias function is separable in ra- 
dial distance r and angular wavenumber £. In this case the 
GI effect can be written in the same way as the tomographic 
shear power spectrum; 



C?/{£) = 



/GI 



drW:-\r)bi{£)bi{r)Pss 



£ 



Sk{r) 



(42) 



This is similar to the (linear/non-linear) alignment model in- 
troduced in Hirata & Seljak (2004) and generalized by Bridle 
& King (2007) and Joachimi & Bridle (2009). The intrinsic 
alignment bias relation, equation (|4ip . assumes that the in- 
trinsic alignment is a non-local convolution of the angular 
density field which depends locally on distance. 

Given we are assuming a model for which we should 
have some motivation to trust, we shall only assume a 
Gaussian-prior on the lA bias functions. In Appendix D we 
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Figure 3. The dark energy Figure of Merit (FoM) for a Euclid-like tomographic weak Icnsing survey marginaUzed over all non-linear 
intrinsic alignment bias functions. We show results as a function of the scatter in radial intrinsic alignment bias, fe/(r) with the angular 
bias fixed at ai,{i)/P = 10"'* (left panel) and angular intrinsic alignment bias, b^ (i) with ai,(r)/P = 10~* (right panel) (effectively 
tracing the along the x and y-axes in Figure [Jjl. The lower (black) lines show the FoM when only the shear signal has a cosmological 
dependence and the GI term is treated as noise. The upper (red) lines show the FoM when both the shear and the radial weak lensing 
dependence in the GI term dependent on cosmology. 
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logio(f^i/i') 

Figure 4. Dark energy Figure of Merit (FoM) contours, for a 
Euclid-like tomographic weak lensing survey marginalized over all 
possible non-linear alignment bias functions, as a function of the 
fractional scatter in radial intrinsic alignment bias (y-axis) and 
angular intrinsic alignment bias (x-axis). Here we have assumed 
the GI term is dependent on cosmology (upper lines in Figure[3ll. 
The FoM reduction is more severe for a given uncertainty in the 
angular intrinsic alignment bias. 



show that the path-integral marginalized covariance for this 
model is given by 

C^i{£) ^ C^A£) + B^.W (43) 

where 

B^,{i) = bj{e) f "draUr)W,f{r)W['J{r)\Pss{e/r;r)f 



olW 



-iGI /«\r^GI, 



+ 1 ^f^ j C^J WCi-rriW. 



(44) 



We assume 6/(r[z]) — C'ip{z)/D{z){l + z) and bi{£) = 1 as 
the fiducial functions for the biases (see Bridle & King, 2007 
and Joachimi & Bridle, 2009 for these definitions) and that 
the scatter in these functions is given by (Ti,(r) and cr(,(£). 
Equation (|18p is again used to derive a Fisher matrix for the 
non- linear alignment bias model (also given in Appendix C). 

Figure |3] shows the dark energy FoM as a function of the 
radial (left plot) and angular scatter (right plot) in the lA 
bias, for our fiducial survey. The lower set of lines assumes 
that the GI terms are independent of cosmology, while the 
upper set of lines includes the cosmological dependence of 
the GI terms. We see that including the GI cosmology de- 
pendence increases the FoM by a factor of around three. The 
angular bias reduces the FoM when ab{(.)/P > 10~^, while 
the radial bias does not affect the FoM until ab{r)/P > 10''. 
This is because the dominant effect of the dark energy equa- 
tion of state is to change the angular part of the lensing 
power spectrum. While the growth rate is also effected by 
changes in the dark energy equation of state, the extra geo- 
metric constraints from the lensing kernel (in both the shear 
and GI terms) help to break any degeneracies introduced. 
As a result the dependence of the FoM on the radial func- 
tional scatter is weak (and even more so when both shear 
and GI contribute to the geometric constraints) . This is also 
in agreement with Joachimi & Bridle (2009) where a large 
but finite-dimensional marginalization has found a similar 
behaviour. 

Figure |4] shows a 2-D contour plot of the dark en- 
ergy Figure of Merit as a function of the fractional scatter 
in radial and angular intrinsic alignment bias. Again it is 
clear that marginalizing over angular, ^-dependent bias has 
a larger impact than the radial bias on the FoM. A con- 
straint in the functional scatter of at/P ~ 10^^, required 
such that the FoM is unaffected, is achievable using current 
techniques applied to future surveys, e.g. Mandelbaum et 
al., 2009; who find a GI functional constraint of ~ 1 using 
WiggleZ. 
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5 CONCLUSIONS 

Path-integral marginalization can account for the impact 
of measuring and marginalizing over unknown functions on 
cosmological parameter estimation. We have constructed a 
likelihood function in cosmological parameter space that 
takes into account marginalization over unknown functions, 
assuming either a flat-prior, were we use information in the 
data itself to measure and remove the effect, or a Gaussian- 
prior, where we use external data to constrain the nuisance 
functions. We assume the distribution in function-space is 
well approximated by a Gaussian functional, but note that 
the non-Gaussian case can be treated with numerical meth- 
ods. 

We use the likelihood to derive an expression for the 
cosmological parameter Fisher information matrix in the 
case that all functions are marginalized over. We have ap- 
plied this to the case of a scale-dependent galaxy bias func- 
tion, relating the galaxy power spectrum to the underlying 
matter power spectrum. As a scale-dependent bias is com- 
pletely degenerate with the matter power spectrum we as- 
sume a Gaussian-prior with variable scatter. We show that 
the marginal errors of Qm, wo and Wa rapidly increase as 
the uncertainty in the fractional power of a scale-dependent 
bias increases above Pb{k)/h^{k) > 10"^, implying we need 
to constrain any scale-dependent bias to a few percent ac- 
curacy. We did not explore redshift-dependent bias models 
here, but the marginalized Fisher matrix implies a similar 
accuracy is needed. In addition we did not explore stochastic 
bias models, which will add a further requirement. There are 
clear ways to constrain scale-dependent bias from redshift- 
space distortions, which also probe the velocity field and 
so are bias-free; by the galaxy bispectrum; or by combining 
galaxy redshift surveys with weak lensing probes of the dark 
matter (e.g. van Waerbeke, 1998; Bernstein, 2009; Matarrese 
et al, 1997; Verde et al., 1998; Guo & Jing, 2009). 

Our second example is in the removal of Intrinsic Align- 
ments effects, in particular the non-local GI effect, on cos- 
mic shear measurements. In addition to nulling (Joachimi & 
Schneider, 2008, 2009) and modelling (Bridle & King, 2007; 
Kitching et al., 2008; Joachimi & Bridle, 2009), path-integral 
marginalization represents a third way to remove and assess 
the impact of this effect. Path-integral marginalization has 
the advantage that it is model-free, and introduces no biases 
in cosmological parameter estimates. For a flat-prior on the 
ellipticity-density relation, the dark energy FoM drops from 
a fiducial FoM = 130, to FoM ^ 100, when we allow the 
GI term to be completely free in both angular and radial 
dependence and measured by the data. 

If we assume the GI ellipticity-density relation is accu- 
rately constrained from external data (e.g. Mandelbaum et 
al., 2009) or complimentary information (e.g. shear-position 
information; Bernstein, 2009) with a Gaussian-prior, and al- 
low the weak lensing effect in the GI term to depend on cos- 
mology, we find that the dark energy constraint is enhanced 
ft-om FoM= 130 to FoM= 300. The FoM is degraded as the 
accuracy drops below a-{P^s)/P~,5 ~ 10~^. The assumption 
of a non-linear alignment model (e.g Hirata & Seljak, 2004; 
Bridle & King, 2007) relaxes this requirement to « 10~^. 

Our results imply the GI effect in weak lensing can be 
removed by path-integral marginalization and that cosmic 
shear can still constrain dark energy to high precision. In 



addition, with no extra assumptions, we have shown that 
the GI effect can itself be used to probe cosmology through 
its dependence on weak lensing. By constraining the GI 
ellipticity-density dependence, using external data, we can 
use the weak lensing dependence in GI to increase our knowl- 
edge of dark energy. 

The path-integral marginalization approach presented 
here is a powerful and general method to assess, control and 
remove the effects of systematics which cannot be simply 
modelled as parametric functions. 
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APPENDIX A: MARGINALIZATION OF SYSTEMATICS IN THE MEAN 

We start with the hkelihood function where the cosmological parameters, 0, and nuisance parameters, ip, both appear in the 
mean, ^{0,tf}). Expanding the nuisance parameter to first order and explicitly marginalizing over variations in the nuisance 
parameters, i/), we find. 



L{D\e) = / dlSil)] 



-l/2(AD-Si,id^^fJ.)C \AD-S^jd^^nf ^-6ll)C^^6ll)*/2 



^(27r)^detC 
First we expand in a Fourier series 



f{x,^i,C) 



cTs I d"fce-5feCfc'^-»fc.(a;- 



^{2^)';AetC^ 



Sid.^l) I ^m^l ^-^k' C sk'* ^-ik' .8 



(45) 



^n^^-lkCk'^-^k.X 



^™^,^-ifc'C.fc" / ^™^g-,{s,/^.fc-fc:)s. 



d"fce-2 



fcCfe' ^(k.aifj,)Cs_ij(djfji.k) -ik.x 



^(2^)"det[C+(a,/i)Q^.(a,/i*)f 
If we carry out the Gaussian integration this reduces to 



1 



(46) 



(47) 



L{D\9) = ((27r)™ det[C + C,'^(9^.m)(9v.,M )])" -^ exp i^--AD[C + C*(a^,/i)(a^^./i)TMD 
We note that this can be encapsulated by a modification of the covariance 

CM = c + ct^id^^fj.){d^^fi'), 

where Cf, is the covariance of the systematic, which acts as an effective prior. This is a generalization of the results from 
Bridle et al., (2002) and Taylor & Kitching (2010) where C = C + agj^fifj,^ . After carrying out a path- integral marginalization 
over the bias functions, the new data covariance matrix is given by 



(48) 



Om — O - 



-*'<*W**')'^^fSf?' 



(49) 



where the derivatives 5fi[tp{x)]/Stpi{x) are functional derivatives. 

This can also be derived using the Woodbury matrix identity (e.g. Press et al., 1990). Using equations (|6]) and pO[) we 
can write the path-integral marginalized likelihood as 



jC = AD C 



C 



Stpi{x' 
where we have remove the constant term, and 



d.'d."^a»M-(.',.").''^'[^(^')] 



HA^') 



C"' AD+ 



i^\-i 



(50) 



(51) 



M.J = F.J + {CD 

where Fij is the path integral Fisher matrix given by equation (|9]). Then using the Woodbury identity the covariance term in 
equation (|50|) can be directly written as equation (|49[l . 

APPENDIX B: NUMERICAL PATH INTEGRATION 

In this paper we make use of the assumption of Gaussianity in function-space. However, in some cases this assumption will be 
a poor approximation to the actual functional distribution. In this case there are two alternative approaches which we review 
here: Numerical Path-Integration and Functional Form-Filling. 

Numerical Path Integration decomposes the functional integration into discrete steps at position Xi, and integrates 
over all possible field values at each position; 



L{D\e) (X 



\{jv[^:^{x)] 

a 

M „ 

n Y[jd^<^i^ 



p[D\e,-^(x)]p[iij{x) 



p[D\e,^{x.)\p[^{x.)]. 



(52) 



Each integral can be evaluated by Monte-Carlo Markov-Chain (MCMC) methods (e.g., Lewis & Bridle 2002). 

Functional Form Filling was developed in Kitching et al., (2009) and uses a Monte-Carlo approach where the bias in 
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the measured parameters is evaluated instead of marginalization. In the Appendices of Kitching et al (2009) it is shown that 
the bias and marginalization procedures are equivalent. Functional Form Filling is also a slight generalization of Monte-Carlo 
path-integration in that arbitrary basis sets are used and the coefficients of the expansion sampled randomly. 

APPENDIX C : PATH-INTEGRAL MARGINALIZATION FOR POWER SPECTRUM INTRINSIC 
ALIGNMENT MODEL 

Starting from Section |4] equation (|34|l we can write the mean signal is the cosmic shear and intrinsic alignment power spectra 
as 

f^^, {£) = a, {£) = Cf^ {() -f Cfl.^M) + Nii , (53) 

where i,j refer to the background source redshifts, Zi, Zj, and we have intrinsic alignments for galaxies in the i**^ redshift bin. 
The covariance of this is 

Gov ([Cg°(£) + C?j!j>dmCFr^{£) + CPri,n.>l{e)]) = 

cf^ie) [cir^w + cS!,y,{i) + cSUw] + [cSf W + cS[,,>„. W] [cft''{£) + cjifi>, w] . (54) 

This includes all covariances between redshift bins pairs. 

Starting from SectionUwe can treat the ellipticity-density cross-power, P^s{i/r, r), itself as an unknown function. Taking 
the functional derivative with respect to the perturbation we find that 

SC^/W ^yOI,. Suit -I') 

In a similar way to Appendix B we can now write the marginalized covariance as 

C^JA^) = C^^i) + [ "dra%Wr,r)W,f{r)W{'J{r), (56) 

Jo 

where ap{£/r,r) = {\SP^/s{£/r,r)\'^) is the functional scatter in P^s and we have assumed the functional covariance is diagonal 
in £ and r. 

By taking the derivative of the mean with respect to a set of cosmological parameters 6a we can now write a general 
expression for the Fisher matrix 

^a5 = i ^ /" ^ {[C^iWr' [9aMM W ^^M. W + ^a^. W ^OMm W]) (57) 



where fj, = {i,j) and u = {l,m) denote pairs of background redshifts. Equations (|60l) and (|57p provide a Fisher matrix for 
weak lensing tomography that include marginalization over all possible GI intrinsic alignment bias functions. 

APPENDIX D: PATH-INTEGRAL MARGINALIZATION FOR THE NON-LINEAR ALIGNMENT 
INTRINSIC ALIGNMENT MODEL 

In the following we will use the following notation (similar to Hirata & Seljak, 2004 and Bridle & King, 2007) 

Cfj'^W^ j " drW?^'='{r)Pss{£/r;r) and C^/ (£) = b' (£) f " drW^f {r)b' (r)Pss{£/r;r) (58) 

Jo Jo 

which can be related to equations (|29[) and (|42|) . We have here introduced two bias functions b'' (r) and b' {£)■ We have assumed 
that the scale and redshift dependent part of a more general bias b^ [r, £) can be separated into two independent functions in 
this case; this will allow us to comment on the relative merit of constraining the bias function in these directions. 

In a similar way to Appendix G let us consider the bias functions as unknown systematics such that we can write them 
as some fiducial function plus some stochastic Gaussian distributed unknown ipi{^) ~ Sb^ (r) and rpj{£) = 3b^ {£). We can write 
down the functional derivatives associated with these perturbations as 

^^rW-5^(,)VK«^(0P.,(fA...) and ^^^'^^^ C?/ {£) Sni£ - £') 

Using equation H58[) and assuming that the cross terms in the nuisance function covariance matrix are zero, {Sbi{r)5b}{£')) — 0, 
and that the auto-covariance are diagonal in r and £ the marginalized covariance becomes 

Cf:U£) = C,A£) + bj{e)J "dra!ir)W[^\r)Wl'J{r) \Psde/r;r)\' + (^) C%\£)C?^{£:) , (60) 

where (l(56/(r)|^) = o\(r) and similarly for £. 
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